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A modified form of the Lam-Bremhorst low-Reynolds-number k-e turbulence model (ref. 1) has 
been developed for predicting transitional boundary layer flows under conditions characteristic of gas 
turbine blades. Since previously reported work (refs. 2,3) has outlined the basic form of the model 
and it's application to zero-pressure gradient flows with free-stream turbulence, this paper will 
primarily describe the application of the model to flows with pressure gradients, and will include 
tests against a number of turbine blade cascade data sets. Also, some additional refinements of the 
model made in recent months will be explained. 


INTRODUCTION 

The difficulty, and yet great importance, of accurately predicting external heat transfer on gas 
turbine blades has stimulated a large amount of research dedicated to understanding and modeling 
transitional boundary layer flows. The primary factors found to influence this phenomenon include 
Reynolds number, free-stream turbulence, pressure gradients and streamline curvature, all of which 
are present on a typical gas turbine blade. Although many approaches can and have been taken in 
modeling this process, our work has focused on exploring and developing the potential of low- 
Reynolds-number (hereafter abbreviated as "LRN") forms of the k-e turbulence model for solving 
this type of problem. 

A variety of different LRN modifications have been proposed in the literature for the purpose of 
extending the validity of two equation turbulence models through the viscous sublayer to the wall 
(see for example reference 4 for a good review). One attractive characteristic of this type of model is 
the seemingly natural process by which boundary layer transition from laminar to turbulent flow is 
simulated without requiring a separate model. However, although some of these models had already 
been used to predict heat transfer on gas turbine blade cascades prior to the initiation of our work 
(refs. 5,6,7), little had been reported documenting predictive capabilities for the less complex case 
of flat plate flow. As a result, our previously reported work began by evaluating two relatively 
popular LRN models with respect to the prediction of transition on flat plates under the influence of 
free-stream turbulence. This work showed that although these models do predict the qualitative 
aspects of transition correctly, the starting location and length over which it occurs is generally not 
well predicted. 

*Work done under NASA Grant NAG-3-579 
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Subsequent to this evaluation a simple modification was proposed which could be applied within 
the framework of a LRN k-e turbulence model, and which was designed to improve the prediction of 
boundary layer transition without affecting the fully turbulent model. The calibration of this model 
and application to flat plate zero pressure gradient flows, also reported earlier, proved to be fairly 
successful. Thus, this work has been pursued and the purpose of this paper is to report on the 
application of the model to predict experimental results of transitional flows with both ffee-stream 
turbulence and pressure gradients, and also to actual turbine blade cascade experiments. 
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NOMENCLATURE AND SYMBOLS 

Empirical parameters in the proposed modification. Correlated as functions of Tu e 
Constants in the two-equation turbulence models. 

Low-Reynolds-number functions. 

Mean and fluctuating static enthalpy 
Total enthalpy 
Turbulent kinetic energy 
Acceleration parameter K^/U 2 (3U/9x) 

Turbulence length scale L=k 1 - 5 /e 
Static pressure 

Production term in the modeled k equation 
Molecular and Turbulent Prandtl numbers 
Reynolds number based on momentum thickness 
Momentum thickness Reynolds numbers at the start of transition 
Turbulence Reynolds number R t = k 2 /Cue) 

Wall Coordinate Turbulence Reynolds number R y = k- 5 y/o 

Stanton number _ _ 

Turbulence intensity, Tu ={ l/3(u' 2 + v' 2 + w' 2 )}- 5 /U 

Mean velocity in the x direction 

Fluctuating velocities in the x, y, z directions 

Stream wise direction coordinate 

Cross-stream coordinate 

Fluid density 

Boundary layer thickness 

Momentum thickness of the boundary layer 

Acceleration parameter A=8 2 (dU/dx)/u 

Dissipation rate 

Molecular viscosity 

Eddy or turbulent viscosity 

Kinematic viscosity, u = |i/p 

Empirical constants in the turbulence models. 


denoting critical 
denoting ffee-stream value 
denoting value at x=0 
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BASIC EQUATIONS AND THE TURBULENCE MODEL 
For the results presented in this work, the following forms of the time averaged continuity, 


momentum, and energy equations were solved. 

a(pu)/a x + a(pv)/a y =o (i> 

pU[3U/3x] + pV[3U/3y] = 9/3y[pau/3y - piTv'] - dP/dx (2) 

pU[9H/ax] + pV[3H/3y] = 3/ay[(p/Pr)aH/ay - phV] + 

a/9y { U[(l-l/Pr)pau/3y - puV] } (3) 

where H=h+U^/2 . Also the fluid is assumed to be a perfect gas with constant specific heat 
A LRN form of the k-e turbulence model is used to determine the turbulent shear stress and turbulent 
heat flux by defining 

-p£V = p t (9U/3y) (4) 

-phV = (p t /Pr t )3h/3y , with Pr t = 0.9 (5) 

and by solving the following two additional transport equations for k and e. 
pU[3k/3x] +pV[3k/3y] = 3/3y[(p + p t /a k )3k/3y]+ p t (3U/3y) 2 - pe (6) 

pU[3e/9x] + pV[3e/ay] = a/a y [(p+p t /a e )3e/ay]+ C el fi(e/k)p t 0U/a y ) 2 (7) 

where |i t = pC^ (k 2 /e) (8) 


The constants C^, C el , C £ 2 , c k , and a e , and the LRN functions fj, f 2 are given in reference 8. The 
function fji has been slightly changed as compared to the Lam Bremhorst form to improve the 
predictions at very low Tu. Previous work had shown that the Lam-Bremhorst model did not predict 
transition for free-stream turbulence levels of lower than about 1 % ( refs. 2,7). Our investigation 
showed that this was due to the particular two parameter correlation for (f^LB=f^.[R t , R y ]) chosen 
by Lam and Bremhorst which, under certain conditions predicted f^»l. Since the function f^ is 
intended to vary only from 0 to 1, this problem could be eliminated by simply restricting the 
magnitude of f^. The details of this are given in reference 8. For the calculations given here, f^ was 
calculated as 


fn. = min( 1 ., .5+.0025*R t , ^,lb) (9) 

The modification introduced to improve boundary layer transition predictions affects the so-called 
production term in the k equation, Pk=M-t(3U/9y) 2 . (Note that the use of the words "production 
term" has been used rather loosely here to refer only to the quantity in the model, not a term in the 
exact form of the k equation.) This is done by limiting the magnitude of Pk before a specified 
stability limit (R.0^125), and then limiting the growth rate of Pk afterward. The details and the 
numerical implementation of this is given in reference 8. However, the basic relationship can be 
expressed as follows; 
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[ aPk/at ] 

max — A*P k + B 


( 10 ) 


The functions A and B have been correlated to the ffee-stream turbulence such that under flat- 
plate conditions the model predicts transition occurring at momentum thickness Reynolds numbers at 
the beginning and end of transition in according with the correlation of Abu-Ghannam and Shaw 
(ref. 9). It should be realized that the behavior shown in Figure 1 is unique to the LRN functions 
used in this model. Somewhat different values must be used with, for example, the Jones Launder 
model (ref. 10) to achieve similar results (see ref. 8). It also reflects somewhat different values then 
reported earlier in reference 2. This is a result of a reevaluation in the high Tu range, and also to the 
minor modification made to the f^ function given in eq. (9). 


COMPUTATIONAL PROCEDURE 

Although details cannot be given here, a brief overview of some of the prac tical aspects of the 
computational procedure is appropriate. Before one can begin any computation, starting profiles and 
boundary conditions must be specified such that the actual problem of interest is solved. For a two 
equation turbulence model, this must include values for k and e in addition to the velocity and 
enthalpy. For all of the results presented here the following procedure was followed. 

Initial Starting Location: The method of Th waites ( see ref. 1 1, pg 3 15 ) was used to compute the 
momentum thickness Reynolds number variation from either the leading edge or the stagnation point 
All calculations were then started at a streamwise location such that the momentum thickness 
Reynolds number was less than about 25. Our previous work showed that this minimized sensitivity 
to variations in initial profiles (ref. 2). 

Velocity; Following the procedure of Rodi and Scheuerer (refs. 6,7), a Pohlhausen polynomial 
representation of the velocity profile was used. This requires an approximation for the local 
boundary layer thickness 8, and an acceleration parameter A=8 2 U'/t). The details of this are given in 
reference 8. The local ffee-stream velocity and velocity gradient was determined from the 
experimental conditions. 

Enthalpy; The total enthalpy in the ffee-stream was assumed to remain constant for all cases 
considered. For flat plate flows the starting enthalpy profile was derived from an approximate 
temperature-velocity profile relationship as per reference 6. For the turbine blade calculations, 
which were started near a stagnation point, this procedure was slightly modified to allow the thermal 
boundary layer Sr, to be different than the velocity boundary layer. This was controlled by 
estimating the stagnation point heat transfer coefficient and relating this to Sr. The details are given in 
reference 8. 

k and e: The ffee-stream values of k and e were found by relating them to the experimentally 
reported values of the free-stream turbulence intensity. For isotropic turbulence, this implies that 

ke=1.5(Tu e U e ) 2 (11) 

The value of e must be found such that the calculated decay in k with distance matches the data when 
the following two ordinary differential equations are solved. 


158 



U e (dk e /dx) = -e e 
U e (d£^dx) = - C e2 (ee 2 /ke) 


(12) 

(13) 


Documentation of these values for each case studied and the comparison with the experimental 
data is provided in reference 8. Since our previous work had shown the value of specifying ke based 
on u’ only, when possible this was always done. 

When the experimental data is not sufficient to determine the dissipation length scale, such as is 
the case with much of the turbine blade cascade data, the calculations assume a small value of e such 
that k remains essentially constant. Note that because of the relationship between k and Tu expressed 
in eq. (10), even if the value of k remains constant, Tu will vary with ffee-stream velocity. 

For convenience, the initial profile specification of k and e is similar to that of Rodi and 
Scheuerer ( ref. 7) except that their parameter "ai" was always assumed equal to unity. Details of this 
are also given in reference 8. 


COMPARISON OF THE COMPUTATIONS TO EXPERIMENTS WITH FAVORABLE 

PRESSURE GRADIENTS 

Taken together, the data sets of Blair (refs. 12,13), and of Rued and Wittig(refs. 14,15,16) 
provide a good variety of well documented experiments where the effects of both acceleration and 
ffee-stream turbulence on transition are represented. Both sets of data also provide the 
experimentally determined variation of Tu over the test sections, a key requirement for accurately 
modeling the transition process. Furthermore, the experiments of Rued have the additional attractive 
characteristic of representing wall to gas temperature ratios that are similar to those experienced on a 
gas turbine blade. 

The experiments of Blair were set up so as to provide a flow with essentially constant 
acceleration, and two such levels of acceleration were explored. The magnitude of the acceleration 
parameter K corresponded to 0.2 x 10 6 , and .75 x 10 A In figure 2, the Stanton number data found 
for the lowest acceleration case is shown for three different turbulence generating grids and 
compared to the computations. Although the Stanton number in the fully turbulent region is under- 
predicted, the location and extent of transition is reproduced fairly well. In Figure 3, results for the 
higher acceleration are shown. Here, the location of transition is reasonably reproduced at the higher 
turbulence level, but deviates somewhat in the lower turbulence case. Compared to the calculation, 
the experiments show the onset of transition to be further downstream, and the extent of the 
transition region to be longer. Furthermore, the extent to which the fully turbulent Stanton number is 
underpredicted has increased. The effect on transition is not unexpected, as we know that 
acceleration has a stabilizing effect on turbulence. However, it is not clear why the fully turbulent 
predictions were low. In our opinion, the most likely possibilities are an inadequate modeling of the 
turbulent heat flux (we used Pr t =.9 = constant), and/or a problem in the dissipation equation in the 
near wall region during acceleration. 

Figure 4 shows the variation over the test section of the acceleration parameter K for three of the 
test cases reported by Rued. Note that these achieve significantly higher values of K than the tests of 
Blair. Figures 5-7 show a comparison of our calculations with the reported Stanton number data for 
these flows at four different levels of free-stream turbulence each. The model is quite successful at 
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predicting the heat transfer for cases 6 and 12, except in the initial region after the unheated starting 
length. One of the important observations made by Rued and Wittig (ref. 13) about their work, was 
that the results did not show a lower limit on Reg s of 163 as is used in the Correlation of Abu- 
Ghannam and Shaw for Tu>7%. This is reflected in these figures by our computations predicting 
the initial transition occurring too late at high Tu. For test case 10, grid two, the model indicates a 
relaminarization, whereas the experiments do not. Also, in grid 1, where transition occurs only 
after the acceleration has stopped, the transition length is underpredicted. In some ways this reminds 
us of the computation of Blair's data shown in figure 3, for grid 2. 


COMPARISON OF THE COMPUTATIONS WITH TURBINE BLADE CASCADE DATA 

Figures 7-10 provide a comparison of the models predictions with the data provided by Daniels 
and Browne (ref. 17), and Hylton et al.(ref. 18). Data at two different Reynolds numbers is shown 
for each blade on both the suction and the pressure side. 

Since the turbulence intensity measurements reported by Daniels and Browne were taken 
somewhat upstream of the test section, figures 7 and 8 show the results of computations assuming 
Tu 0 =3.5% and 3 % for comparison purposes. As can be seen, the location and extent of transition, 
and the resulting heat transfer is very well predicted for these cases. The only significant variation 
between the data and the computation occurs at the higher Reynolds number in regions downstream 
of transition. On the suction side Starting at S/L= 4-.5, the flow experiences a region of adverse 
pressure gradient Previous research has documented the failure of 2-eq. LRN models such as that 
of Lam-Bremhorst to correctly calculate the near wall turbulence length scale in adverse pressure 
gradient flows, resulting in an over prediction of the skin friction and heat transfer (see ref. 19). 

For comparison, a computation is shown where the dissipation equation was modified in line with a 
suggestion of Launder (ref. 20) in the following manner. 

C el =1.44*max( 1, L/L^) 


where 

L=k 15 /e > L max =2.7y ( 14) 

The results of this additional modification show an improved prediction of the heat transfer on the 
suction side without influencing the transition predictions. 

In figures 9 and 10, comparisons with the data of Hylton et al. are shown. These calculation 
show similar trends as was pointed out on the Daniels and Browne data, i.e., an excellent prediction 
of the lower Reynolds number data, but an over prediction of the heat transfer in the turbulent 
regions on both the suction and pressure sides for the higher Reynolds number case. Also, the 
modification given in eq. (14 ) is shown to reduce the predictions on the suction side down to match 
the data. However, there is a clear overshoot in the heat transfer predictions for run 145, which 
cannot be explained with reference to the dissipation equation correction. The reason for this 
difference may be related to the effects of curvature in die recovery region after the release of 
curvature. In the C3X blade the radius of curvature is small until about S/arc - 2 . 
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CONCLUSIONS 


Tests of the transition model proposed against a wide range of flows with pressure gradient have 
shown that even without additional modifications, the model is quite satisfactory for predicting much 
of the data. However, it appears that the computations might be improved in the future by 
considering the following; 

(1) As the relationship between Re© and Tu at high Tu is clarified experimentally, improvements 
could be incorporated into the correlations for A and B, and possibly Reec- 

(2) The length over which transition occurs during favorable pressure gradient conditions might be 
more accurately modeled by decreasing the value of A and/or B in some appropriate manner. 

(3) Improvements in the LRN functions relative to the fully turbulent flows with pressure gradients. 

(4) Appropriate incorporation of curvature effects into the transition modifications. 
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FIGURE 1 Variation in "A" and "B" 
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RE 2 Comparison of the calculations with the lower acceleration test data of Blair 



FIGURE 3 Comparison of the calculations with the higher acceleration test data of Blair 
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FIGURE 6 Comparison of the calculations with the data from case 10 of Rued 




Grid 4, T 14= 1 0% 


□ □ D □ □□ 
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S/L (suction side) 

FIGURE 8 Comparison of the calculations with the data from Daniels and Browne 
on the suction side of the blade L=50.44 mm 
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S/L (pressure side) 

FIGURE 9 Comparison of the calculations with the data from Daniels and Browne 
on the pressure side of the blade L=50.44 mm 
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Tu =6.5% 


Calc without adverse PG corr. 


u m m m m a ■ 




°nn Oj 


• Run 145, Re2=2.49E+6 
n Run 149, Re2=1.51E+6 


S/arc (suction side) 


FIGURE 10 Comparison of the calculations with data on the suction side from 
runs 145 and 149 of the C3X blade of reference 18. arc=.18 m 
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Run 145, Re2=2.49E+6 
Run 149, Re2=1 .51E+6 


S/arc (pressure side) 

FIGURE 1 1 Comparison of the calculations with data on the pressure side from 
runs 145 and 149 of the C3X blade of reference 18. arc=.14 m 
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